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Abstract 

qj ■ The charge asymmetry of the g, h, and k Dalitz plot parameters for K — > -k 7r°7r° 



decays has been measured with 35 GeV/c hadron beams at the 70 GeV IHEP 
accelerator. The g, h, and k values obtained appear to be identical for K decays 
within the errors quoted. In particular, the charge asymmetry A g = {g + — 9~)/(g + + 
g~) of the slope g is equal to (0.2 ± 1.9) • 10~ 3 . 

X 
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1 INTRODUCTION 



Observation of direct CP violation in neutral kaon decays [1,2,3] gives argu- 
ments to search for a similar effect in charged kaon decays. For example, this 
effect can manifest itself as a charge asymmetry of the Dalitz plot parame- 
ters for K^ — > 7r ± 7r°7r° decays. These parameters are coefficients in a series 
expansion of the squared module of the matrix element [4]: 



\M(u,v)\ 2 oc 1 + gu + hu 2 + kv 2 , 



1 Corresponding author. E-mail address: kozelov@mx.ihep.su 
Preprint submitted to Eur. Phys. J. C 7 February 2008 
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Fig. 1. Experimental layout: M - magnets, Q - quadrupoles, CM - tuning mag- 
nets, K - collimators, S - scintillation counters, C,D - threshold and differential 
Cherenkov counters, BH - beam hodoscopes, AC - anticoincidence counter, H - 
scintillation hodoscopes, GEPARD - electromagnetic calorimeter. 

where u and v are invariant variables. 



Theoretical estimates of the charge asymmetry of the Dalitz plot slope g for 
K ± — ► 7r ± 7r°7T decays are uncertain and range from 10~ 6 to 10~ 3 [5,6,7,8]. 
In the majority of the experiments, only g + or g~ was measured [4,9]. From 
these studies it follows that Ag = g + — g~ = 0.066 ± 0.017. It is very un- 
likely to expect direct CP violation at this level, and one can assume that the 
above mentioned difference is due to the underestimation of the systematic 
uncertainties. 



K — »■ 37r decays have been studied simultaneously for both K + and K mesons 
in [10,11,12]. Ford et al. [10] found A g = -0.0070 ±0.0053 for K± -»■ 7T ± 7T ± 7r q: 
decays. In the experiment [11], A g = 0.0019 ±0.0123 was measured for K^ — ► 
^i^o^o ,-[ eca y S# Preliminary analysis of our experimental data [12] based on 
a fraction of statistics yielded A g = —0.0003 with a statistical error of 0.0025 
and a systematical uncertainty below 0.0015. In this paper we report our final 
results on the charge asymmetry of the Dalitz plot parameter measurements. 



2 EXPERIMENTAL SETUP 



The experiment was carried out with the TNF-IHEP facility [13] (Fig. 1) at 
the 70 GeV IHEP accelerator. Unseparated 35 GeV/c positive and negative 
hadron beams used for kaon decay studies are produced by 70 GeV protons 
in the external 30 cm Al target. Scintillation counters S1-S4 and beam ho- 
doscopes BH1-BH4 are used to monitor beam intensity and to measure parti- 
cle trajectories and beam profiles. The typical particle flux was 4xl0 6 per 1.7 
second spill. 
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Fig. 2. Efficiency of the differential Cherenkov counters vs argon pressure. 

Kaons are selected with three threshold C1-C3 and two differential Dl, D2 
gas Cherenkov counters (Fig. 1). The admixture of unwanted particles under 
the kaon peak was substantially below 1% (Fig. 2). The threshold counters 
are also used to select 10 GeV/c electrons to calibrate GEPARD calorimeter. 

About 20% of kaons decay in the 58.5 m long vacuum pipe located downstream 
of the BH4 hodoscope. The flanges of the vacuum pipe have thin Mylar win- 
dows in the path of beam particles. The 3.6 m diameter exit flange is made of 
4 mm thick (0.23 Xq) stainless steel. The probability of a high-energy photon 
to convert into an e + e~ pair in this flange is equal to 0.16. 

Kaons which pass through the decay pipe are detected by the anticoincidence 
counter AC. The BH5 beam hodoscope placed behind the calorimeter is used 
for a high precision measurement of the beam position at the setup end. The 
BH5 hodoscope operates in the counting mode and hence detects all beam 
particles. 

The products of kaon decays are detected by three scintillation hodoscopes 
H1-H3 [14] and the GEPARD electromagnetic calorimeter. Each hodoscope is 
made of two X,Y octagonal planes with 3.85 m distance between the oppo- 
site octagonal sides (Fig. 3). The plane is divided into half-planes with 256 
elements each. The cross section of the hodoscope elements is 14 x 12 mm 2 
and their length varies from 0.7 to 1.8 m. Scintillation light is detected by 
FEU-84-3 photomultiplier tubes. 
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Fig. 3. General view of the scintillation hodoscope H. 




Fig. 4. The GEPARD calorimeter. 
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Fig. 5. M 77 distribution for if — > -k it decays. 

The GEPARD is a sampling lead-scintillator calorimeter. It contains 1968 
cells with 76 x 76 mm 2 cross section (Fig. 4). Each cell consists of 40 alter- 
nating layers of 3 mm Pb and 5 mm scintillator. The total radiation length 
is 21 X . Scintillation light is collected onto FEU-84-3 photomultiplier tubes 
using wavelength shifting light guides. The GEPARD calorimeter was cali- 
brated by irradiating each cell with 10 GeV/c electrons at the beginning of 
data taking and by using K ± — > ir^it reconstructed events collected during 
the experiment. Both methods yielded consistent results. The 7r° mass resolu- 
tion is equal to 12.3 MeV/c 2 (Fig. 5). 

The Level 1 trigger is formed according to the logic formula 



Tl = SI ■ S2 ■ S3 ■ S4 ■ (Dl + D2) ■ CI ■ C2 ■ C3 ■ AC. 



The Level 2 trigger uses information about energy deposition in the GEPARD 
calorimeter [15]. For this purpose the calorimeter is divided into 16 trigger ele- 
ments. The Level 2 trigger is formed if the energy deposition exceeds 0.8 GeV 
in at least three trigger channels. 



The stability of the beam and detector parameters was carefully monitored 
during the data collection. To reduce the systematical uncertainty in the mea- 
surement of the charge asymmetry of the Dalitz plot parameters the beam 
polarity was reversed every day. Figs. 6-7 illustrate the stability of the exper- 
imental setup operation. 
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Fig. 6. Time variations of the probability P n to reconstruct n tracks (P n is averaged 
over ~10 5 events). 
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Fig. 7. Time variation of the mean energy per event detected in the GEPARD 
calorimeter (Ec a i is averaged over ~10 5 events). 

3 EVENT RECONSTRUCTION AND SELECTION CRITERIA 



The K± 



Tr^TT 



event selection starts by finding energy clusters in the 



GEPARD calorimeter. The coordinates of the cluster centers and the X and Y 
coordinates measured by H1-H3 hodoscopes are used in track reconstruction. 
To reduce the combinatorial background, only tracks with three or four hits 
in each X and Y projection are selected. Then the vertex position of the K^ 
decay is calculated using the reconstructed tracks. A track is considered to be 
associated with a kaon decay if the hypothesis of its intersection with beam 
axis has a confidence level of 5% or more and the decay vertex position is 
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Fig. 8. x 2 distribution for K — » it 7r°7r° decays (histogram - simulation, circles - 
experiment). 

inside the fiducial volume of the decay pipe. In addition selected events have 
to satisfy one of the following criteria: 

• five clusters with energies above 1 GeV are found and each track is associ- 
ated with one of these clusters; 

• four clusters with energies above 1 GeV are found and one of the tracks is 
not associated with these clusters. 



These criteria are applied because there is a substantial probability for the 
gamma from 7r° decays to convert into e + e~ pair in the exit flange of the decay 
pipe (see Section 2), and charged pion energy deposition in the calorimeter 
could exceed the threshold value of 1 GeV. 

Events passing this preliminary selection are subjected to a kinematic fit that 
allows one to resolve ambiguities due to the combinatorial background (for 
example, to associate one of the tracks with the charged pion) and to calculate 
the Dalitz plot variables. Altogether 21 measured parameters are used in the 
fitting procedure: the energies and the coordinates of four clusters associated 
with gammas, the kaon mean energy and the parameters of the kaon and 
pion tracks. The parameters of the clusters are corrected for the transverse 
profile of the electromagnetic shower and for the spatial nonuniformity of the 
calorimeter. The energy of the charged pion is the only unknown parameter. 

Seven constraints are imposed on the fitted parameters: four equations of the 
energy-momentum conservation, two equations for the effective masses of the 
gamma pairs and a required intersection of kaon and charged pion trajectories. 
The decay vertex coordinates are not fixed. The parameters are found by the 
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Fig. 9. Number of events vs kinematic variables for K 
- simulation, circles - experiment). 



ir it decays (histogram 



minimization of the functional with constraints using the method of uncertain 
Lagrange multipliers and iteration technique. The iterations are stopped when 
the relative changes of all fitted parameters at the last iteration are less than 
10~ 5 . For each event all possibilities to associate one of the tracks with charged 
pion and 7 pairs with 7r°'s are considered. The combination with the least x 2 
is used. Fig. 8 shows the \ 2 distributions for the data and simulated events. 
Events with x 2 > 20 are rejected since in this region the data exceeds the 
number of the simulated events due to the high background level. Simulation 
shows that this x 2 cut reduces statistics of the wanted events by 28% only but 
decreases background level considerably. 

The experimental setup operation was simulated using a Monte Carlo (MC) 
method with the GEANT 3.21 code. The setup geometry is described in de- 
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Fig. 11. it, v resolutions (u,v - "true" variables, u',v' - measured variables). 

tail, and the data obtained in the experiment were taken into account. Among 
these data are the calibration coefficients for each channel of the calorime- 
ter, the dependence of the hodoscope efficiency on the particle coordinates 
and correlations between kaon spatial and angular coordinates and its mo- 
mentum. Fig. 9 presents different distributions of the experimental data and 
simulated events for the K ± — ► tt^tt decays. From this figure it follows that 
there is a good agreement between experimental and simulated distributions. 
Fig. 10 shows the acceptance of the setup, and Fig. 11 demonstrates the u,v 
resolutions averaged over the Dalitz plot. 

The x 2 probability P(x 2 ) for data and simulation is shown in Fig. 12. Events 
with P(x 2 ) > 0.1 are selected for further analysis since in this region there 
is a good agreement between the experimental and simulated data. To check 
the Level 2 trigger conditions the energies corresponding to each of the trigger 
channels are calculated. The event is accepted if the number of the channels 
with energy above 1 GeV is greater than two. This cut rejects only a few 
K^ — > 7i ± n iT events (see the last row of Table 1), but it is important for 
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adjust the simulation code, and to estimate the systematic uncertainties. 
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The final data sample is comprised of iV + = 278398 and N~ = 341015 events. 
Table 1 shows the fraction of events rejected by each cut and the cumulative 
efficiency. 

Table 1 



Selection criteria 


Fraction (%) of events 


rejected by the 
cut 


passed this and all 
previous cuts 


> 1 track is reconstructed in H1-H3 


4.4 


95.6 


Position of the decay vertex is inside 
the fiducial length of the decay pipe 


31.2 


65.7 


Number of clusters and tracks 
corresponds to the K — > tt tt 7t° decay 


93.2 


4.48 


X 2 <20 


82.2 


0.80 


P{ X 2 ) > 0.1 


26.4 


0.59 


Level 2 trigger is ok 


0.2 


0.59 



After all selection cuts for K^ — > 7r ± 7r°7r° decays, an admixture of background 
remains. The background sources are other modes of kaon decays, interactions 
of beam particles in the material along the beam line and overlapping of events 
due to the finite time resolution of the detectors. Simulations of these processes 
demonstrate that the main contribution to the background comes from the 



10 



1.5 

1 

0.5 



-0.5 

-1 

-1.5 



n 


n 


n 


□ 


□ 


□ n 








\c 


l_l uuuuuu 

□□□anoo 
□nnnnna 


LJ 


□ 


□ - 

On 


n 














II 




II II 


— 














- - [ - - -|- - - 1 - - - 


□ 














1 




1 l| ! 
















1 




I ' I 




m 




_ 














i i | 






— 










- - J- - - 




II 






□ 


~ , 










1 




1 


















ll 




U 




















II 




in 





















--]!- 




3E 


-- 








= 


- 






— 


1 


= 


Ho 


o 




mo 








O 


oo 


Oa 
□ ■ 


□ □ 

□ n 


□ □ 


□ □an 
□□:□-□ 

□ □:□□ 

□ □:□ n 

□ □ a 

II 1 1 1 II 


no 
□ n 


l_IO 




■ 




□ 




□ 


□ 




ii 



0.5 



1.5 

1 

0.5 



-0.5 

-1 

-1.5 



-n n n n n nn a 

cyan ocratrcna 



K 




□ □□:□□ no □□ 

□ .a □:□ □.□:□.□ □ 

□ □ □:□ □ □:□ □ □ 

□ □ □:□ □ □:□ n □ 
o □ □ ;□ □ □ | □ ■ 



□ no 

□ □.a 

□ □ ■ 



1.5 2 2.5 0.5 1 1.5 2 2.5 

\V'\ \V'\ 

Fig. 13. Distribution of events in the Dalitz plots for K^ — > 7r ± 7r 7r decays. 



K ± -»■ n^n (0.21% ) and X ± -> ^ir+iT (0.03%) decays. This contribution 
does not depend on the sign of the kaon charge and hence does not cause a 
false charge asymmetry of the Dalitz plots. The background level from other 
sources is less than 0.01%. 



The finite energy resolution of the calorimeter results in a noticeable (~ 10%) 
probability of the wrong combinations of 7's reconstructed into 7r°'s and the 
hodoscope inefficiency can cause a reconstruction of a false track (~ 5%). Both 
these effects are taken into account in the event simulation and it is found 
that their influence on the charge asymmetry of the Dalitz plot parameters is 
negligible. 



4 RESULTS 



4-1 Difference of the Dalitz plot parameters 



The difference of the Dalitz plot parameters for K^ decays is estimated by 
minimizing the following functional form: 



X \Ag, Ah, Ak) = E K ' ~ l ~ aiA9 ~/ lAh ~ llAk? : 



(2) 



hj 



v 



%/AH 



<7% 



y,Z 



-4r + -=■ I , where n,-,- is the number of events in the 



' l 3 n-/ M - ' "« ' l 3 \n+ ' nil' ""^ "'"V 

13 IN \ 13 <-3 / 

i-th, j'-th Dalitz plot bin with u', v' measured coordinates (Fig. 13), and ay, 
Pij, and 7jj are coefficients defined in (A. 2) (see Appendix) and calculated by 
MC. The values of Ag, Ah, and Ak as well as the elements of the correlation 
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matrix are given in (3): 



Ag = -0.0009 ±0.0067, 
Ah = -0.0007 ±0.0062, 
AA; = -0.0014 ±0.0017, 



^1.00 0.93 0.35^ 



V 



1.00 0.32 
1.00 



(3) 



The errors shown are statistical only. The x 2 /ndf is 319/(279-3) = 1.16. 



Fig. 14 shows the rj(u') 



j^n-./N- aIld r J'( U ') 






ratios of the nor- 



malized event distributions in Dalitz plots for K — ► n n n decays projected 
on the u' and \v'\ axes. The Tij(u') ratio versus u' for different \v'\ intervals are 
presented in Fig. 15. 

Since some theoretical models predict that CP violation in K ± — ► 3tt decays 
can be associated with the charge asymmetry of the slope g only, Ag is also 
estimated assuming Ah = Ak = which is in agreement with (3). With this 
assumption we find: 



Ag = 0.0002 ± 0.0024 



X 2 /ndf = 319/(279 



1.15. 



(4) 



The g, h, and k parameters appear to be equal for kaons of different signs 
within statistical uncertainties. But this does not guarantee the identity of 
the event distributions in the corresponding Dalitz plots. In order to check 
the identity of u', \v'\, and («', \v'\) distributions independently of the ma- 
trix element form (1), the Kolmogorov nonparametric criterion was used. 
This analysis provided the following results: the probabilities that u', \v'\ and 
two-dimensional («', \v'\) distributions for K + and K~ are indistinguishable 
are equal to 32.4%, 85.4% and 55.2% correspondingly. The two-dimensional 
(«', \v'\) distributions were compared using the modified Kolmogorov criterion 
from the HBOOK code which works with histograms. 



4-2 Systematic errors 



All measures are taken to assure that K + and K~ beams have identical pa- 
rameters. Nevertheless the average angles Ax and Ay of beam particles with 
respect to the nominal beam axis and the mean kaon energies in the positive 
and negative beams could differ by AAx = 5 firad, AAy = 7 firad and AE 
= 50 MeV. Simulations show that these uncertainties result in the following 
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systematic errors: 

S A (Ag) = 0.0004, 5 A (Ah) = 0.0003, 5 A (Ak) = 0.0001, 
8 E (Ag) = 0.0006, S E (Ah) = 0.0004, 5 E (Ak) = 0.0001. 

Two methods were used to evaluate systematic uncertainties of Ag, Ah, and 
Ak connected with g, h, k errors (see (2) and (A.1),(A.2) in Appendix): MC 
and the analytical method assuming ideal resolutions in u and v : u = u', 
v = v' (see (A. 3), (A. 4)). Both methods gave the same results: 



S(Ag) 

Ag 
5 (Ah) 

Ah 
S(Ak) 



^(0.2 • 5g) 2 + (0.6 • 5h) 2 + (1.6 • 5k) 2 , 



(0.5 • 5g) 2 + (1.0 • 5h) 2 + (1.5 • 8k) 2 , 
(0.4 • Sg) 2 + (0.5 • 5h) 2 + (2.9 • 5k) 2 . 



Ak 

Using experimental data [4] one can obtain: 

5(Ag) 5(Ah) 5(Ak) 

-V^ = °- 014 > -^TT 1 = 0.024, -A— ■+■ = 0.019. 

A^ Ah Ak 

Other possible sources of the systematic errors include the time variations 
of the calorimeter calibration coefficients and the hodoscope efficiency, the 
influence of the Earth magnetic field on the particle beams of different polarity, 
the difference in the 7r + and ^interactions with matter, and the difference 
in composition and intensity of the positive and negative beams. The total 
contribution of these factors to the systematic errors does not exceed 1 • 10 -4 . 

We also investigated that varying the minimum energy of 7's, the minimum 
and maximum energies of charged pion, the value of the x 2 confidence level 
and the number of reconstructed tracks did not appreciably change the results. 
The results remain also unaffected if the bins located at the boundary of the 
Dalitz plot are not used. Ag, Ah, and Ak can also be found by minimizing the 
functional for the differences of the Dalitz plots. The obtained results agree 
with (3) and (4). 

Thus, the final estimates of the systematic errors are 

5(Ag) = 7-10~ 4 , 5{Ah) = 5 • 10~ 4 , 5(Ak) = 1.4 • 10~ 4 . (5) 

The systematic errors are approximately an order of magnitude less than the 
statistical errors given in (3). 
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CONCLUSIONS 

The differences Ag, Ah, Ak of the Dalitz plot parameters have been measured 
for the K^ — > ir^ir^n® decays using the TNF-IHEP facility. The studies were 
performed in the 35 GeV/c positive and negative hadron beams at the 70 GeV 
IHEP accelerator. Frequent changes of the beam polarity allow one to minimize 
the systematic uncertainties of the experiment. Our results show that the event 
distributions in the Dalitz plots for K + and K~ decays are indistinguishable 
and that the Ag, Ah, and Ak are equal to zero within the errors quoted 
(3), (5). Assuming Ah = Ak = we find: 

Ag = 0.0002 ± 0.0024(stat.) ± 0. 0007 (syst.). 

To find A g = Ag/(g + + g~) it is assumed g + = g~ = 0.652 [4]: 
A g = 0.0002 ± 0.0018(siot) ± 0.0005 (syst.). 

This is the most accurate estimate of the charge asymmetry of the Dalitz plot 
slope for the K ± — > 7r ± 7r 7r° decay. 
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APPENDIX 



A Calculation of the differences of the Dalitz plot parameters 



In accordance with (1), the probability density function can be expressed in 
the form: 



f(u',v> 



J G • (1 + gu + hu 2 + kv 2 )dudv 



D 



J J G ■ (1 + gu + hu 2 + kv 2 )dudvdu'dv' ' 

DD 



where u, v and u',v' are the "true" and measured Dalitz variables, and G = 
G(u, v, u', v') is a function that depends on the efficiencies of the detectors and 
data processing. Integration is performed over the kinematic boundary of the 
Dalitz plot. This relation can be rewritten as: 

., , A a + gb + hc + kd 



1 + gu + hu 2 + kv 2 



where a = a(u', v') — - J G dudv, b = b(u', v') — - J u ■ G dudv, c = c(u', v' 

£ d £ D 

1 fu 2 

D D DD 



ju 2 -Gdudv, d = d(u',v') = -jv 2 -Gdudv, u = - J J u ■ G dudvdu'dv', 
u 



2 = - J J u 2 ■ G dudvdu'dv' ', and v 2 = - J J v 2 ■ G dudvdu'dv', and e = J J G dudvdu'dv'. 

£ DD £ D D DD 

The u, u 2 , and v 2 are the mean values of the Dalitz variables and their squares, 
and e is the total efficiency of the experiment (including the event reconstruc- 
tion efficiency) for |M| 2 =1. 

Let's introduce the following notations: 

g = {g + + g~)/2, h=(h+ + h~)/2, k = {k + + r)/2, 

Ag = g + -g~, Ah = h + - h~ , Ak = k + - AT. 



f-(u',v') 

Ah, and AA; and neglecting their quadratic terms, we obtain 



Expanding the f -\^, ,v v ,\ ratio of the normalized Dalitz plots into series in Ag, 



f + (u' v') 
r(u', v') = J - / ; « 1 + a(u', v')Ag + (5{u' , v')Ah + 7 ( u ', v')Ak, (A.l) 
f-{u',v') 



where 
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a(u', v') = [b — au + h(bu 2 — cu) + k{bv 2 — du)]/D(vf, v'), 
P{u', v') = [c-a^ + g[cu - brf) + k{crf - ch?)]/D(u', v'), (A.2) 

-y(u', v') — [d — av 2 + g{du — bv 2 ) + h(du 2 — cv 2 )]/D(u', v'), 

D{u\v') = (1 + g u + hu 2 + kv 2 ) ■ (a + gb + hc + kd). 

In the case of u' = u and v' = v, ("ideal" resolution) formulae (A.1),(A.2) 
have the form: 



= i + Aq(u, v)Ag_ + Bp(u, v)Ah + C (u, v)Ak 

(1 + g u + h u 2 + k v 2 ) ■ (1 + g u + h u 2 + k v 2 ) 



where 



A (u, v ) = u — u + h(uu 2 — u 2 u) + k{uv 2 — v 2 u), 
B (u,v) = u 2 — u 2 +g(u 2 u — uu 2 ) + k(u 2 v 2 — v 2 u 2 ), (A. 4) 

Cq(u, v ) = v 2 — v 2 + g{v 2 u — uv 2 ) + h{v 2 u 2 — u 2 v 2 ). 
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